Magnetism of Ru and Rh thin films on Ag(OOl) substrate 
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In a very recent x-ray magnetic circular dichroism experiment concerning with Ru and Rh im- 
purities and metal films on Ag(001) substrate, no local magnetic moments were displayed in direct 
contradiction with previous theoretical works. It is thought that there can be three main reasons 
for this inconsistency: relaxation, alloying and many-body effects. Some of the above-mentioned 
systems are studied by using a first-principles method in which relaxation and alloying are taken 
into account, even so magnetism is still obtained. For low-coverage systems, high magnetic moments 
in both Ru (~ 2.49 /is) and Rh (~ 2.00 /zb) are obtained. Naturally, as the coverage is increased 
the magnetic moments are approached to zero. Also, it is noticed that the relaxation distances 
are increased by magnetism, which in turn is decreased by alloying. The behavior of the magnetic 
properties is explained in terms of Stoner model. 
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I. INTRODUCTION 



Low-dimensional systems involving transition 
metals (TM) is the subject of an extraordinary 
research activity which is driven by both funda- 
mental and technological interest. The complex- 
ity of this field is illustrated by considering that 

frequently the experimental^'2'^£&£&2 ' 10 i 1 1 and 
theoreticali^ii^^^^o^^ works are in 

contradiction each other. For example, the theoretical 
possibility of ferromagnetism in Ad TM thin films 
on noble substrates was originally reported by Zhu^ 
and Eriksson^ then several experimental works were 
performed"^: some of those with the intention of 
creating^ that kind of systems in the laboratory and 
others with the idea of searching for magnetism 4,5 . 
However, magnetism was not found experimentally. 
This lack of accord was regarded as a consequence of 
that neither relaxation nor pseudomorphism were taken 
into account in the theoretical models. In this sense, 
Wu and Freeman^ 5 -, using a full potential method which 
took into consideration relaxation by means of total 
energy minimization, presented a comparison between 
magnetic and paramagnetic states as a function of the 
Ru(Rh)-Ag interlayer distance. In the same direction 
Bliigeli^ ignoring the force minimization reported the 
dependence of magnetism as a function of the layer 
coverage. Afterward, Tureki 7 - et al. using an ab initio 
method and making allowance for mixing with the Ag 
atoms, obtained the magnetic moment as a function 
of the coverage and the mixing with the Ag substrate. 
Besides, experimental works done by Chang^ et al. 
and Bcckmann 9 et al. presented the possibility of 
clusters formation in the surface, and due to this the 
vanishing of magnetism. Following these experimental 
works Stepanyuki^. et al. studied small and mixed Ad 
clusters on Ag(001) surface, and in according with other 
theoretical works, they found magnetic moments in Ru 
and Rh atoms. Much more recently, Honolkai et al., 



using x-ray magnetic circular dichroism obtained no 
magnetic signal in Ru and Rh thin films on Ag(001) 
in contradiction with theoretical work o 13 i 14 ' 16 . Also, 
Chadoii et al. using STM images reported absence of 
ferromagnetism in a study of Rh atoms deposited on 
Au(lll) substrate. 

The purpose of this paper is to give an advance in 
the comprehension of the subtle interplay between mag- 
netism and geometry in Ru and Rh thin films on Ag(001). 
This specific problem is expected to be particularly in- 
teresting from a fundamental point of view due to the 
current inconsistency between theory and experiment^. 
In particular, it seems that two out of three mean rea- 
sons for this disagreement have been ruled out. 

The remainder of the paper is organized as follows. In 
the next section a brief outline of the ab initio method 
used for the calculations is given. The Results obtained 
for paramagnetic state and the magnetism in Ru and 
Rh thin films on Ag(001) fee substrate are presented and 
discussed in section 3. There have been considered effects 
that in previous works have not been reckoned. Finally, 
section 4 summarizes the main conclusions. 



II. THEORETICAL METHOD 

We have used the ab initio method called Plane Waves 
Self Consistent Field2&. This method, based on den- 
sity functional theory^ 7 ., is restricted here to use the 
generalized gradient approximation of Perdew-Burkc- 
Ernzerho^ (GGA-PBE) functional. 

The studied systems are modelled using slab geome- 
try. Our model consists of various geometrical configu- 
rations divided in supported (noted by O) (Fig. 1 a.l), 
alloyed (Fig. 1 a. 2) and buried (Fig. 1 a. 3) thin films 
on Ag(001) substrate. Supported atoms are modelled by 
R^Ei-^, where x = 0.25 and 0.50 and R(E) represents 
the Ru or Rh atoms (empty space). Alloyed atoms are 
modelled by R x Agi_ x , where x is the same as above. 
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FIG. 1: Scheme of the geometrical configurations, (a) The 
systems are modelled in 3 cases, (a.l) supported, (a.2) alloyed 
and (a.3) buried atoms. Positions in the unit cell consist of 4 
sites labelled by a, b, c and d inside the dashed line, see (b) 
and (c). The case (b) is for R0.25E0.75 where only the site a 
is occupied by the R atom (Ru or Rh) the others sites being 
empty space or Ag atoms. The case (c) is for R0.50E0.50- The 
R atoms are located at sites a and b whereas c and d are 
empty sites or Ag atoms. 
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FIG. 2: Paramagnetic local density of states (LDOS) for Ru 
(a) and Rh (b) atoms (only d component) in Ru and Rh thin 
films on Ag(001) substrate. We present the LDOS for Ru and 
Rh atoms for each system in alloyed and supported systems 
with x = 0.25, x = 0.50, 1 and 2 ML (surface (Surf.) and 
Interface(I)). Fermi level is dotted line. 



of Honolkai et al. 



A. Paramagnetic state 



The unit cell has 4 inequivalent atoms per plane. These 
geometrical configurations are showed in Fig. 1(6) and 
1(c). Also we studied 1 and 2 monolayers (ML). The geo- 
metrical configurations consist of 7 metallic layers, where 
5 layers correspond to Ag(OOl) substrate and 1 adlaycr 
on each side of the slab. These adlayers represent the 
supported, alloyed atoms and full monolayers depending 
of the case. The metallic layers are separated by 7 layers 
of empty space, which is sufficent to prevent the interac- 
tion between the slabs and to vanish the charge at the 
central layer of empty space. All atoms are located at 
the ideal positions of a (001) slab with lattice parameter 
of Ag(OOl) (4.09A), however we determined the atomic 
position using force minimization. The calculation are 
performed using the Monkhorst-pack scheme to define 
the k points for each slab with a grid of 16 x 16 x 1. A 
cut-off energy of 36 Ryd (489eV) was used for the plane 
waves expansion of the pseudowave function (560Ryd for 
the charge density and potential). The intcrlayer dis- 
tances are relaxed until the absolute force is less than 
0.001 Ryd/a.u. 



III. RESULTS 

The aim of this paper is determinate the electronic and 
structural properties of Ru and Rh thin films on Ag(001) 
using an ab initio method, the studied systems are pre- 
sented in Fig. [TJ and 1(2) full monolayer(s), in order to 
explain the absence of magnetism in experimental results 



The Stoner model^L is able to reveal the magnetic or- 
der in 3d transition metals and their alloys, with few 
values of the system to consider. The needed values to 
indicate the magnetic order are the paramagnetic Local 
Density of States (LDOS) and the Stoner parameter (is). 
Therefore, the Stoner model is able to describe ferromag- 
netism, weak ferromagnetism, supcrparamagnetism and 
paramagnetism. 

We report the paramagnetic LDOS (only d component) 
for Ru and Rh atoms, see Fig. O 

Following the Stoner model we have the next equation: 

S = IsNnm(Ef) 

where Is is the Stoner parameter and Nnm(Ef) is the 
value of paramagnetic LDOS at the Fermi level. The 
I s values are 0.2992 and 0.3264 (in cV) for Ru and Rh, 



TABLE I: Values of S factor for each system. The value for 
atoms supported (O), alloyed (A), buried (B) surface (Surf) 
and interface (Inter). The letters a b and c are the atomic 
positions. 





Ru 


Rh S/A/B 


0.25 a O/A 


2.91/1.96 


1.46/0.47 


0.50 a O/A/B 


3.05/2.18/1.28 


2.23/1.27/0.54 


1.00 O/B 


1.08/1.08 


1.43/1.12 


2.00 Sur (Inter) 


0.64(0.67) 


0.48(0.55) 
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TABLE II: Difference of energy (in meV/Atom) between mag- 
netic and non magnetic state for Ru and Rh thin films on 
Ag(OOl), in the magnetic and non magnetic state. The cap- 
ital letter represent Supported (O), Alloyed (A) and Buried 
(B). 



Ru 


Rh 


0.25 O/A -29.86/-12.71/ 
0.50 O/A/B -69.12/-32.20/-11.15 
1.00 O/B -33.67/-20.71 
2.00 O 0.00 


-49.97/-25.82/ 
-121.05/61.23/-13.85 
-86.06/-61.24 
105.72 


TABLE III: Magnetic moments (in hb) for Ru and Rh atoms 
in Rh and Ru thin films on Ag(001) substrate. The capital 
letter represent supported (O), alloyed (A) and buried (B). 


Ru 


Rh 


0.25 O/A/S 2.49/1.86/ 
0.50 O/A/B 2.56/2.12/1.48 
1.00 O/B 2.07/1.65 
2.00 O 0.01/0.00 


2.00/1.80/0.85 
1.99/1.79/1.15 
1.57/1.40 
1.75/1.51 



TABLE IV: Distances (in A) between R atoms and Ag(001) 
substrate, in the magnetic (zm) and non magnetic (znm) 
case. The O and A represent the Supported and Alloyed, re- 
spectively. 





Ru 

O 

ZNM / ZM 


A 

ZNM / ZM 


Rh 

O 


A 


0.25 


1.60/1.72 


1.91/1.94 


1.62/1.85 


1.91/2.03 


0.50 


1.73/1.80 


1.90/1.96 


1.74/1.91 


1.91/2.05 


1.00 


1.98/1.99 




1.67/1.74 





respectively^. We report the obtained S values in Table 

in 

For low coverage, the LDOS presents a narrow and 
big peak below the Fermi level, this peak decreases when 
the coverage increases. The LDOS presents a big value 
at Fermi level which satisfies the Stoner criterion for low 
coverage, and when the coverage increases the Stoner cri- 
terion is not satisfied. The S factor decreases when the 
coverage increases. 

We can see in Table fl] that the S factor (and conse- 
quently the magnetic order) changes as a function of the 
coverage and the type of system (O, A, or B). It means 
that the system has a magnetic phase transition from 
strong ferromagnetism (S >> 1) to weak ferromagnetic 
(2 > S > 1) to superparamagnetic (S f» 1) and finally 
paramagnetic (S < 1). 

We can see in Table U that the S factor for Rh is al- 
ways less than 2, except for the case of supported surface 
ordered alloy, where it is slightly bigger than 2. It means 
that for the case of Rh thin films on Ag(001) present 
ferromagnetism only for the supported surface ordered 
alloy. In other cases the system presents non strong fer- 
romagnetism. 

In the case of Ru we can numerate as following: the 



S factor is bigger than 3 for x = 0.50 for Supported 
atoms; 5=2.91, 2.18 and 1.96 for x = 0.25 Supported, 
x = 0.50 Alloyed and x = 0.25 Alloyed atoms, re- 
spectively; 5=1.28, 1.08 and 1.08 for x = 0.25 Buried, 
x = 1.00 Supported and Buried ML, respectively. Ru 
systems have a phase transition with differences respect 
to Rh systems. 

In both cases (Rh and Ru) for 1 ML, the S factor is 
slightly bigger than 1, it means that 1 ML satisfies the 
Stoner criterion for superparamagnetism or weak ferro- 
magnetism, and 2ML do not satisfy the Stoner criterion 
(S<1). 

We can conclude that alloying and buring (interdiffu- 
sion) reduce the S factor and consequently change the 
magnetic order. We may find ferromagnetism for very 
low coverage and for supported atoms, and when the cov- 
erage increases the magnetic order changes. 



B. Magnetic state 

The possible reasons for vanishing magnetism are al- 
loying, relaxation and many body effects, in this section 
we report the propicrties of the magnetic state, in all 
cases the magnetic state is the ground state. First we 
define the difference of magnetic energy as the difference 
between the magnetic and paramagnetic state, it is given 
by 

AEj\i ag = Em — Enm 

where Em and E^m are the energies in the magnetic and 
non-magnetic state respectively. The values of AEMag 
are presented in Table HT1 We report in Table HVl the dis- 
tance (in A) between R atom and the Ag(001) substrate 
in the magnetic (zm) and non magnetic (znm) case. Also 
we present the magnetic moment for Ru and Rh atoms 
in Table ULTl for each system. 

In almost all cases the magnetic energies are bigger 
than the thermal energy (25 meV). The alloying reduces 
the magnetic moment, this is due to the neighborhood of 
Ag atoms. The interdiffusion or buring reduces the mag- 
netic moments. The magnetism increases the distance 
to the Ag surface, in agreement with magneto-bulk cou- 
pling. 

The Stoner criterion is not satisfied when S < 1, how- 
ever we find magnetic moments. One possible explana- 
tion is the enhanced susceptibility that polarize the Ru 
and Rh atoms. 

We cannot comparate our results with previous the- 
oretical works because the distances used here are ob- 
tained by fully relaxation or force minimization. 



C. Cluster formation 

Experimental works have presented the cluster 
formation^ and interdiffusion^.^. In order to take into 
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TABLE V: Energy differences (in meV/atom) respect to the 
ground state (noted by 0.00) in the magnetic and non mag- 
netic state. 



System 


NM 


Mag 


Ru 






Surface Ordered Alioy 


230.20 


189.10 


Surface 2D Aggregate 


93.89 


82.74 


Buried Ordered Alloy 


114.69 


93.20 


Buried 2D Aggregate 


10.35 


0.00 


Rh 






Surface Ordered Alloy 


114.74 


65.76 


Surface 2D Aggregate 


117.73 


66.53 


Buried Ordered Alloy 


37.91 


24.70 


Buried 2D Aggregate 


30.63 


0.00 



account these phenomena, we cosider the energy for 2D 
surface aggregate and 2D buried aggregate by the follow- 
ing equation: 

E 2 d = ^{Er 1ai0 + E Agi00 ) (1) 

where Er 100 and EA gi 00 are the energies for 1ML of R 
and Ag on Ag(001). Using Eq. |T]) we determine the dif- 
ferences of total energy for each system, and reported in 
Table IVl In both cases Ru and Rh the system more stable 
is the magnetic buried aggregate, and in both cases the 
S factor is almost 1, it means that this type of aggregate 
satisfies the Stoner criterion for supcrparamagnctism or 
weak ferromagnetism, in other words the Stoner criterion 
is not satisfied for strong ferromagnetism (S >> 1). 



IV. CONCLUSION 

Following the experimental work of Honolka 1 about 
Ru and Rh overlayers films on Ag(001) substrate^, we 
have performed a theoretical study of the Ru and Rh 
thin films on Ag(001), by means of an ab initio method. 
In order to explain the experimental results, i.e. the ab- 
sence of magnetism, we considered two out of the three 
mean possible causes: relaxation and the alloying effects. 
However, magnetism occurred. It seems then that the 
origin of the discrepancy between the mentioned experi- 
ment and current theory must be sought in other causes, 
for example, many-body effects, spin fluctuations, steps 
among others. The alloying reduces the magnetic en- 
ergy by a factor of 2, and the magnetic moments per 
atom as well. The magnetism increases the distance be- 
tween the atoms, in agreement with the magneto-bulk 
coupling. We have used the Stoner model to reveal the 
magnetic order and we have found that the magnetic or- 
der is a coverage function, and that the systems have a 
magnetic phase transition in coverage function. For very 
low coverage the system is ferromagnetic, and when the 
coverage increases the magnetic order goes to weak ferro- 
magnetic, superparamagnetic and finally paramagnetic, 
successively. By considering the surface aggregate forma- 
tion and buried aggregate, we determine the difference of 
total energy and the ground state, in both cases is the 
magnetic 2D buried aggregate, and this aggregate does 
not satisfy the Stoner criterion. The two main conclusion 
arc: (1) the Ru and Rh thin films have a magnetic phase 
transition as a function of the coverage and (2) the 2D 
buried aggregate is a more stable system. 
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